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We calculate the form factors and the coupling constant in the pD* D* vertex in the framework 
of QCD sum rules. We evaluate the three point correlation functions of the vertex considering both 
p and D* mesons off-shell. The form factors obtained are very different but give the same coupling 
constant: gpD'D* ~ 6.6 ± 0.3f . This number is 50 % larger than what we would expect from SU(4) 
- - I estimates. 
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' Charmonium production is a very useful source of information in heavy ion collisions. The knowledge of the J/ip 
production rate can improve our understanding of these collisions and help us to know if there was a "color glass 
condensate" in the initial state. Charmonium production is very sensitive to the existence and to the properties of 
, the intermediate "quark gluon plasma" All the interesting effects happening in these initial and intermediate 
^Li phases can be blurred by interactions in the final stage of these collisions, when charmonium states interact with 
^ other comovers such as pions, p mesons and nucieons, which form a hot hadronic gas. Since these interactions occur 
^ at an energy of the order of magnitude of the temperature, (~ 100 — 150 MeV), their study has to be made with 
, non-perturbative methods. These can be QCD sum rules 12J, quark models and the effective Lagrangian approach 
' 0i 3 ■ This last approach has been developed for almost ten years now and a great progress in the understanding of 
^ ' the interactions of charmed mesons with light mesons and nucieons has been achieved. Part of this progress is due 
^ I to a persistent stud y o f the vertices involving charmed mesons, namely D*Dn [5, S], DDp^, DDJ/i) H], D*DJ/il) 
00 ; 0*0*1: [IE 112, D*D*J/^ [13, DsD*K, D*DK [I3 and DDuj More specifically, it is very important 

, to know the precise functional form of the form factors in these vertices and even to know how this form changes 
00 ' when one or the other (or both) mesons are off-shell. This careful determination of the charm form factors has been 
'""J ' done bit by bit over the last seven years in the framework of QCD sum rules, which are the best tool to give a first 

' principles answer to this problem. 
^-H ' Understanding charmonium production in heavy ion collisions would be already a good reason to study of hadronic 
charm form factors. However, since 2003, due the precise measurements of B decays performed by BELLE, BES and 
BABAR, this subject gained a new relevance. In B decays new particles have been observed, such as the Z3s,/(2317) 
and the X(3872). These particles very often decay into an intermediate two body state, which then undegocs final 
k>( state interactions, with the exchange of one or more virtual mesons. As an example of specific situation where a precise 
rS , knowledge of the pD* D* form factor is required, we may consider the decay A'(3872) J/ip + p. As suggested 
' in [13, this decay proceeds in two steps. First the X decays into a. D - D* intermediate state and then these two 
particles exchange a D* producing the final J/ip and p. This is shown in Fig. lb and If of [13- In order to compute 
the effect of these interactions in the final decay rate we need the pD*D* form factor. 

In the present paper we calculate this form factor with QCDSR. The pD*D* vertex is similar to the J /ipD*D* 
vertex treated in [l3|. As before, because there are three vector particles involved, the number of Lorentz structures 
is very large and we have to choose a reliable one to perform the calculations. Here we introduce the pole-continuum 
analysis and impose the pole dominance as a criterion to reduce the freedom in the choice of the Borel parameter. 
In the next section, for completeness we describe the QCDSR technique and in section III we present the results and 
compare them with results obtained in other works. 

II. THE SUM RULE FOR THE pD* D* VERTEX 

Following our previous works and especially Ref. [13, we write the three-point function associated with the pD*D* 
vertex, which is given by 
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for an ofF-shell p"*" meson, and: 

for an ofF-shcU D*~ meson. The general expression for the vertices ([1]) and ^ has fourteen independent Lorentz 
structures. We can write each r^Q,^ in terms of the invariant amphtudes associated with each one of these structures 
in the foUowing form: 

(p^ , p'^ , q'^)gf^aP^ + Te {p^ , p'^ , q'^)guap'f, + T7 {p^ ,p'^,q^ )PtJ.PuPa + Ts (p^ ,p^,q'^ )p'f,PuPa 

+'^9{p^,p'^,q'^)Pp.plPa +Tio{p^,p''^,q'^)pf,p^p'^ + Tii{j)'^ ,p''^ ,q^)p'^p'^pc + Ti2{j)'^ ,p'^ ,q^)p'^pup'a 

+^l3(/,p'^ + ^l4(p^p'^ (3) 

Equations ((1]) and ^ can be calculated in two diferent ways: using quark degrees of freedom -the theoretical or 
QCD side- or using hadronic degrees of freedom -the phenomenological side. In the QCD side the correlators are 
evaluated using the Wilson operator product expansion (OPE). The OPE incorporates the effects of the QCD vacuum 
through an infinite series of condensates of increasing dimension. On the other hand, the representation in terms of 
hadronic degrees of freedom is responsible for the introduction of the form factors, decay constants and masses. Both 
representations are matched invoking the quark-hadron global duality. 



A. The OPE side 



In the OPE or theoretical side each meson interpolating field appearing in Eqs. ((T|) and ([2]) can be written in terms 
of the quark field operators in the following form: 

jP^ {x) ^ d{x)-^yu{x) (4) 

and 

{^) = c{x)l^Ax) (5) 

where m, d and c are the up, down and charm quark field respectively. Each one of these currents has the same 
quantum numbers of the associated meson. 

For each one of the invariant amplitudes appearing in Eq.([3]), we can write a double dispersion relation over the 
virtualities p^ and p'^, holding — —q^ fixed: 

^,b^p'^Q^) = -4 r rf^r z = i,...,i4 (6) 

^ Js^,„ A„.„ {S~p^){u-p'^) 

where pi{s,u,Q^) equals the double discontinuity of the amplitude Ti{p^ ,p'^ ,Q^), calculated using the Cutkosky's 
rules. The invariant amplitudes receive contributions from all terms in the OPE. The first one of those contributions 
comes from the perturbative term and it is represented in Fig. [TJ 



D* 




D'- c ]j,o p+ u jj,- 



FIG. 1: Perturbative diagrams for the p off-shell (left) and D* off-shell (right) correlators. 



We can work with any structure appearing in Eq. ([3]) , but we must choose those which have less ambiguities in the 
QCD sum rules approach, which means, less influence from the higher dimension condensates and a better stability 
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as a function of the Borel mass. We have chosen the g^aQi/ structure. In this structure the quark condensate (the 
condensate of lower dimension) contributes in the case of D* meson off-sheh. 
The corresponding perturbative spectral densities which enter in Eq. ([6]) are 



p('')(s,u,Q2) 



27r\/A 



^ (2m2)j {——\+2{J-I)^-{2Tal-u-s)-D 



for p off-shell, and 



U—S — t\ ( B — A\ ^, TT , „ 9x ^ 



(7) 



(8) 



for _D* off-shell. Here A = A(s, m, t) = + + - 2st - 2su - 2tu, s ^ , u ^ p'^, t = -Q^ and A, B, D, I and J 
are functions of (s, i, u), given by the following expressions: 



cos 9 



B = 27ri^cos( 



7r|fc| 



1 - cos^ 



cos — ko 
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D ^ -n\k\ (l-cos^ 



1 - cos 9] cos( 



where 



Po 
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|fc| 



s + u — t 

ir2 2 
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cos 6* 



fcn = 



4s' 2yi 
u + rjml — 2pgfco . 
2b-' 1^ 



(9) 

(10) 
(11) 

(12) 
(13) 



with r\—\ for p off-shell and yy = — 1 for D* off-shell. 

The contribution of the quark condensate which survives after the double Borel transform is represented in Fig. [2] 
for the D* off-shell case and is given by 



(14) 



where (qq) is the light quark condensate. For the p off-shell there is no quark condensate contribuition. 

We expect the perturbative contribution to dominate the OPE, because we are dealing with heavy quarks. For this 
reason, we do not include the gluon and quark-gluon condensates in the present work. 



D* 




p (uu) D* 

FIG. 2: Contribution of the uu condensate to the D* off-shell correlator. 
The resulting vertex functions in the QCD side for the structure g 1^0(1)1^ are written as 



{s — p'^)(u ^ p' ) 



(15) 
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for p ofF-shell and 

47r-^ Jo {s-p^){u-p' ) 

for D* off-shell, where, as usual, we have already transferred the continuum contribution from the hadronic side to 
the QCD side, through the introduction of the continuum thresholds sq and uo [13] ■ 



B. The phenomenological side 



The pD*D* vertex can be studied with hadronic degress of freedom. The corresponding three-point functions, 
Eqs. ^ and ([2]), will be written in terms of hadron masses, decay constants and form factors. This is the so called 
phenomenological side of the sum rule and it is based on the interactions at the hadronic level, which are described 
here by the following effective Lagrangian Q 



C 



(17) 



from where one can extract the matrix element associated with the pD* D* vertex. The meson decay constants, fp 
and /d* , are defined by the following matrix elements: 



and 



\D*{p))=mD*fn^e'g,ip). 



(18) 



(19) 



where and e^, are the polarization vectors of the p and D* mesons respectively. Saturating Eqs. H]) and ^ with 
the p and two D* states and using Eqs. and we arrive at 



rip) 



(P2 +rnl.){Q^ + ml){P'' + m]j,) 



{p + pT 5^ " + (V - pY g" - {2p - p') 



(20) 



when the p is off-shell, with a similar expression for the D* off-shell: 



-gu^' 



pD'D' 



P'uP'y' 



-9pp' 



m 



D' 



gaa' 



(p + pTg 



{2p-pr g"^ 



{2p'-pYg<-' 



(21) 



The contractions of v' and a' in the above equation will lead to the fourteen Lorentz structures appearing in 
Eq. ([2]). We can see from Eq. (|20p that the form factor g^p]j, o*iQ^) same for all the structures and thus can be 

extracted from sum rules written for any of these structures. The resulting phenomenological invariant amplitudes 
associated with the structure gap,((l)v are 



i%\p'\Q' 



)-g%,D-AQ') 



V2fl,fpml,,mp{2 



-) 



(F2 + ml,)iQ^ + m2)(P'2 + mj,,) 



for the p off-shell, and 



pD'D' 



fh,fpml,mp{Q^ + Aml,) 
' (P2 + m2,.)(Q2 + ml,){P'^ + ml)V2ml, 



(22) 



(23) 



for D* off-shell. 
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In order to improve the matching between the two sides of the sum rules we perform a double Borel transformation 
[l3| in the variables = —p^ and P'^ = — p'^ M'^, on both invariant amplitudes F and Tph- Equating the 

results we get the final expressions for the sum rules which allow us to obtain the form factors 3p5*_D* iQ'^) appearing 
in Eqs. (|22p - (|23p . where T is p or D* . In this work we use the following relations between the Borel masses IVP and 
M'2: = for a D* off-shell and AP = M'^ for a p off-shell. 



mc(GeV) 


rriD' ( GeV) 


mp(GeV) 


/i5-(GeV) 


/p(GeV) 


{qq){Geyf 


1.35 


2.01 


0.778 


0.240 


0.161 


(-0.23)^ 



TABLE I: Parameters used in the calculation. 



III. RESULTS AND DISCUSSION 



Table U shows the values of the parameters used in the present calculation. We used the experimental value for 
fp = with 5p = 4.79 [l^ , and took fo' from ref. [l^. The continumn thresholds are given by sa — {m + A^)^ 
and uq = {m + A„)^, where m is the mass of the incoming meson. Using ~ A„ ~ 0.5 GeV for the contimmm 
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FIG. 3: gio'D'iQ'^ = 1-0 GeV^) as a function of the Borel mass M^. 



thresholds and fixing = IGeV^, we found a good stability of the sum rule for .gp'^*^)* for Af^ in the interval 
1 < < 9 GeV^, as can be seen in Fig. [3] Within this interval wc need to choose the best value of the Borel mass 
to extract the coupling constant of the vertex. It is well known in QCDSR that if we choose a too small value of the 
Borel variable M^, then the sum rule will be dominated by the pole, but the convergence of the OPE is poor. On the 
other hand, if is too large, then the OPE convergence is good but the sum rule is dominated by the continuum. 
The best value of the Borel mass is the one with which both criteria are reasonably satisfied. 

In Fig. [4] we show the pole contribution (solid line) and the continuum contribution (dashed line) divided by their 
sum as a function of Borel mass. In the case of p off-shell, we see that the pole contribution is bigger than the 
continuum one in the Borel window 1 < < 3 GeV^. The best Borel mass seems to be JVP ~ 2.0 GeV^. Also for 
this value of the coupling as is aproximately 0.2 and this suggests that perturbative corrections are small. 

In the case of gp^Jjj, the interval for stability is 1 < < 10 GeV^, as can be seen in Fig.O In order to choose the 
Borel mass we proceed in the same way as before and we analyse the pole and continuum contributions. As indicated 
in Figiniin the window 0.5 < M"^ < 1.5 GeV^ the pole contribution dominates. We choose = 1.5 GeV^. 
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M'(GeV') 



FIG. 4: Pole (solid line) and continuum (dashed line) contribuition to <?pfl*£,*(Q^ = lGeV^,M^), as a function of the Borel 
mass M"^ . 
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FIG. 5: g%'}o, (Q^ = 1 GeV^) as a function of the Borel mass. We show the perturbative contribution (dashed line), quark 
condensate contribution (dotted line) and total (solid line). 

Having determined we calculated the dependence of the form factors. We present the results in Fig. [3 
where the circles correspond to the g^po* o'iQ'^) form factor in the interval where the sum rule is valid. The squares 
are the result of the sum rule for the 5p2*£).((5^) form factor. 



In the case of an off-shell p meson, our numerical results can be fitted by the following exponential parametrization 
(shown by the dotted hue in Fig. [7]): 

4^^.^.(Q2) = 5.22e-QV"o. (24) 
As in our previous works 0, H, H, define the coupling constant as the value of the form factor at 



7 




0,2- 



0,0 0,5 1,0 1,5 2,0 2,5 3,0 3,5 4,0 4,5 

M'(GeV') 



FIG. 6: Pole versus continuum contributions to (?po.o.(Q^ ~ 1 GeV^) as a function of the Borel mass M^. 
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FIG. 7: gp'o*£)» (circles) and Qp^J^t (squares) QCDSR form factors as a function of Q^. The solid and dashed lines correspond 
to the exponential parametrizations of the QCDSR data with the two forms mentioned in the text. 

— — m^, where mm is the mass of the ofF-shell meson. Therefore, usmg = ^"^p in Eq the resulting 

coupling constant is: 

QpD'D' = 6.55 . (25) 

For an ofF-shell D* meson, our sum rule results can also be fitted by an exponential parametrization, which is 
represented by the dashed line in Fig. [71 

5S*.l..(Q^)=4.95e-«Vi3.33. (26) 

Using = — m|,. in Eq (pS)) we get: 

gpD'D' = 6.70, (27) 
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in a good agreement with the result of Eg. ((25)) . 

In order to study the dependence of our results with the continuum threshold, we vary As „ between 0.4 GeV < 
As.u < 0.6 GeV in the parametrization corresponding to the case of an off-shell p. As can be seen in Fig. [U this 
procedure gives us an uncertainty interval of 6.40 < gpO'D* < 6.92 for the coupling constant. 




2 3 4 

tf(GeV') 



FIG. 8: Dependence of the form factor on the continuum threshold for the p off-shell case. The dashed curve corresponds to 
As,i, = 0.5 GeV, the solid one to As,u = 0.6 GeV and the dotted curve to As,„ = 0.4 GeV. 

Concluding, the two cases considered here, off-shell p or £)*, give compatible results for the coupling constant, 
evaluated using the QCDSR approach. Considering the uncertainties in the continuum thresholds we obtain: 

gpD-D* =6.6 ±0.31 . (28) 

This generic value of the coupling constant can be easily related to the couplings of the specific charge states. From 
Eq. P?|) we arrive at: 

9p-D'0D'+ gp+D*OD'+ .„„^ 

gpD'D- — — — = — — = -gpOD*+D'+ = gpOD'OD'o 

From Eqs. ([M]) and we can also extract the cut-off parameter, A, associated with the form factors. We get 
A ~ 1.64 GeV for an off-shell p meson and A ~ 3.65 GeV for an off-shell D* . The cut-off values obtained here follow 
the same trend as observed in Refs 0, [1, EO] '■ the value of the cut-off is directly associated with the mass of the 
off-shell meson probing the vertex. The form factor is harder if the off-shell meson is heavier. 

As for the value of gpD" d* , this coupling has not been discussed in the literature as much as those involving the 
J/f/; and there are only few works presenting estimates for it. The starting point in these estimates is always the 
SU(4) symmetry. According to SU(4) we should expect: 

gj/ipD'D' ~ gj/ipDD (30) 
gpD'D' — gpDD (31) 



and 



gpD'D- - "^gj i^D'D' (32) 



From our previous works [1, [T^ we find that in QCDSR Eq. ([50]) is satisfied. However, from [7] and from the 
present work we conclude that Eq. (pij) is not satisfied, since gpw d* ~ 6.6 ± 0.26 whereas gpOD = 3.04 ± 0.30 (20| . 
Eq. ([32)1 is not true either because gjj^D'D' = 6.2 ± 0.9. These relations are violated at the level of 50 %. This is 
not surprising since the mass difference starts to play an important role when we go from the heavier vector mesons 
to p. As for the absolute value, the existing estimates, used in [1, lead to gpWD- — 2.52. Our result is a factor 
two larger. 
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